Intensity limits for stationary and interacting multi-soliton complexes 
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We obtain an accurate estimate for the peak intensities of multi-soliton complexes for a Kerr-type 
nonlinearity in the (1+1) - dimension problem. Using exact analytical solutions of the integrable 
set of nonlinear Schrodinger equations, we establish a rigorous relationship between the eigenval- 
ues of incoherently-coupled fundamental solitons and the range of admissible intensities. A clear 
geometrical interpretation of this effect is given. 



Pulse or beam evolution under the action of compet- 
ing linear dispersion and nonlinear self-action is a fun- 
damental phenomenon which can be observed in optical 
fibers, nonlinear crystals or gases, atomic Bose-Einstein 
condensates, etc. The corresponding wave properties can 
often be described by the nonlinear Schrodinger equa- 
tion (NLSE). In the (l-l-l)-dimensional case, if the ho- 
mogeneous medium has a Kerr-type nonlinear response, 
this equation is integrable and its solutions can be 
found analytically. They describe a nonlinear superpo- 
sition of localized waves, called solitons, and radiation 
modes. Such analysis allows us to reveal some general 
features of nonlinear localized waves. 

Localized solutions of the NLSE are among the main 
objects of interest for telecommunication and informa- 
tion processing applications. All possible localized so- 
lutions of integrable systems are composed of solitons. 
Hence it is important to understand their interactions 
and collective properties. For a single NLSE describing 
the evolution of a coherent field, soliton interactions are 
phase-sensitive, so that constructive or destructive non- 
linear interference can be observed. Since the profiles of 
interacting solitons are distorted according to the nonlin- 
ear superposition principle, large variations of the peak 
intensities can occur. Indeed, it has been demonstrated 
that the peak intensity can vary by a factor of N"^ for 
N interacting solitons |^. In some sense, phase sensi- 
tivity can be "amplified" due to nonlinearity, so that it 
is greater than in the case of linear interference between 
mutually coherent sources. 

Interactions of fundamental solitons can also be phase- 
insensitive if they correspond to incoherently coupled 
components. Then, fundamental solitons can be com- 
posed into a stationary multi-soliton complex (MSG). To 
have a rough idea of an MSG, we can use the follow- 
ing analogy. If we think of solitons as elementary par- 
ticles, then we can imagine an MSG as a more compli- 
cated structure, such as an atom. An optical incoherent 
soliton is an example of an MSG. These objects have 
been studied extensively both experimentally and theo- 
retically during recent years (see Refs. |^^||] and refer- 
ences therein); some other examples of MSGs are listed 



in Ref. 0|. The question of the maximal intensity of 
stationary and interacting MSGs is of primary interest 
for applications. However, to the best of our knowledge, 
this problem has not been addressed yet. In this work, 
we derive a rigorous relationship between the parameters 
of incoherently-coupled solitons and the range of admis- 
sible intensities for MSGs. 

The evolution of spatial incoherent solitons in a 
medium with a Kerr-type self-focusing nonlinearity can 
be described by a set of NLSEs for mutually-incoherent 
components [p|: 
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where Um is the normalized amplitude of the m-th com- 
ponent (1 < m < M), M is the number of components, 

/ = [""jl^ t^*^ total intensity, z is the coordinate 

along the direction of propagation, and x is the trans- 
verse coordinate (we are considering a (l-l-l)-dimensional 
case). 

Eqs. (|l|) are completely integrable by means of the in- 
verse scattering technique (1ST) |^ and, therefore, any 
localized solution can be represented as a nonlinear su- 
perposition of solitary waves. Every fundamental soli- 
ton (labeled j) is characterized by a complex eigenvalue 
kj = rj -\-ifJ.j-, and a shift in the coordinate plane {xj, Zj). 
The simplest single-soliton solution can be written as: 
p^Tj ^ where (3j = rj (xj — /ij Zj ) , 

■Zj) 



-\- (r^ — ^'^)zj/2, and (xj,Zj) 



(x Xj 



Uj{x,z) — rj sech(/3j) 

are the shifted coordinates. The solution for bright 
MSGs consisting of M incoherently-interacting funda- 
mental solitons can be found by solving a set of linear 
equations |g|Jlot]: 
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where Cj = Xj exp (/3j -I- ijj) and the coefhcients Xj de- 
fine the relative coordinate centers for the fundamental 
solitons. The seemingly straightforward process of in- 
verting the set of linear equations results in a complicated 
expression. 
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Using the technique described in Ref . , the solution 
of Eqs. (|^) can be obtained in a compact exphcit form. 
However, it is not possible to find a general analytical 
expression for the maximum intensity, and it can only be 
determined by numerically solving transcendental equa- 
tions. Therefore, a different approach is needed to make 
an analytical estimate for the peak intensity levels. In or- 
der to do this, we first multiply Eqs. (H) by u*/{kj — i/i), 
add the complex conjugate, and sum over the fundamen- 
tal soliton numbers 1 < j < M (we note that such a 
transformation was suggested earlier as an interme- 
diate step for deriving bright and dark soliton solutions 
on a background). Finally, we obtain the following in- 
equality 



M 

E 



< 1, 



(3) 



which is valid for any (real) ^. This key result makes it 
possible to estimate the limitations on the total inten- 
sity. In particular, we conclude that the peak intensity is 
limited from above, i.e. Imax{z) = maxa;/(x, z) < linf, 
where the upper boundary is 



inf 



c[r2 + (^, -^)2], 



(4) 



and /X is chosen to minimize the value of linf- We will 
now consider several examples illustrating the physical 
meaning of Eqs. (^) and (||). 
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FIG. 1. (a) Geometrical illustration of the maximum in- 
tensity criterion for a stationary MSG; black dots mark the 
eigenvalues of fundamental solitons in the (r, jj,) parameter 
space, and the radius of the shaded half-circle determines the 
maximum amplitude, (b) The intensity profiles corresponding 
to upper plots; shading shows allowed ranges of peak intensi- 
ties. 

Stationary MSC. Let us first analyse the properties of 
a single MSC, composed of several solitons with iden- 
tical velocities, In such a case, the optimal choice 
of the free parameter in Eq. (^ is = /i^, and we 
have linf — maxj(r|). Note that this value can be 
interpreted as the minimal squared radius of a circle 
which has its center at the point (0, ^) in the parameter 



space (r, n) and which contains all the soliton eigenval- 
ues within it [see Fig. |l|(a)]. Thus linf is proportional 
to the area {-kR^ /2) of this semi-circle. On the other 
hand, we note that for a stationary MSC the ampli- 
tude profiles satisfy a self-consistent eigenvalue problem. 



-.U„ 



21 Urn = 0, where the functions 



Um — Um'3xp{—i^m) are real. Then, since solutions 
should be localized in the transverse direction (x), we 
conclude that Imax > hup = max^ (r|/2) . This means 
that the variations of the peak intensity are strictly lim- 
ited by the largest eigenvalue in an MSC, as illustrated 
in Fig. |l| Note that Imax — hnf if there is only one 
fundamental soliton (Fig. |^, left). The peak intensity 
decreases if several solitons compose an MSC, but al- 
ways remains above the lower limit (Fig. 0, right). This 
happens despite the fact that the individual fundamental 
soliton intensities are always superimposed, i.e. destruc- 
tive interference is not possible, in contrast to the case 
of coherent interactions. The observed decrease of total 
intensity underlines the complicated nature of the non- 
linear superposition phenomenon and occurs because the 
profiles of individual solitons are strongly distorted due 
to the nonlinear self-action effect. 




FIG. 2. (a) Geometrical illustration of the maximum in- 
tensity criterion for interacting MSCs: black dots mark the 
eigenvalues of fundamental solitons in the (r, fj) parameter 
space, and radius of the shaded half-circle determines the 
maximum amplitude, (b) Dependencies of the peak inten- 
sities, corresponding to the upper plots, on the propagation 
distance; dashed line shows the maximum possible value. 



Interacting MSCs. Our general results can also be ap- 
plied to the case when solitons have different velocities, 
/i„. In other words, we can estimate the peak intensity 
changes during the collision of several MSCs. As follows 
from the form of Eq. (^, the limiting value /„/ depends 
only on the maximum eigenvalue in each of the colliding 
MSCs. Again, Eq. (||) has a clear geometrical interpre- 
tation: the optimal value of ji in Eq. (^) must be chosen 
to minimize the area of a half-circle which has its cen- 
ter at the point (0, /x) and which contains all the soliton 
eigenvalues. Two examples, corresponding to a collision 
between a single fundamental soliton and an MSC with 
different velocities {^n), are shown in Fig. ||. When the 
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relative velocity is small (Fig. g, left), the intensity at 
the impact area of the collision decreases. We have al- 
ready observed this effect for a stationary MSG when 
the relative velocity is zero. However, for larger relative 
velocities, the peak intensity increases (Fig. 0, right). 

In order to understand the differences in the interac- 
tion pattern, we study the dependence of the peak in- 
tensities on the relative velocities of the colliding MSCs. 
To illustrate the key features, we consider interactions 
of two identical MSCs with a relative velocity 2/zi [see 
Fig. ^(a)]. Although Eq. (^) can be used to obtain an 
estimate for the peak intensities, as in the previous ex- 
amples, we find that the results are not optimal for large 
fii. To obtain a more accurate estimate, we recall that 
Eq. (^) is satisfied for all (real) fj, simultaneously. We 
now choose /i = ±/zi, add up the corresponding in- 
equalities together, and obtain the following upper limit, 
Itnf = maxj [r|(r| -I- 4^i)/(r| -I- 2^1)]. We thus see that 
linf ifJ'i +Qo) 2maxj(r|), which is a simple sum 
of the upper bounds for individual MSCs. This result 
means that interaction of solitons with large relative ve- 
locities is weak, and their intensities are added together 
similarly to the linear case. This conclusion is confirmed 
by numerical simulations, the results of which are pre- 
sented in Fig. ||(b). We also illustrate the evolution of 
the intensity profiles at small and large relative velocities 
in Figs, ^(c) and (d), respectively 
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FIG. 3. (a) Location of the eigenvalues of fundamental 
solitons in the (r, /x) parameter space (black dots) . (b) Ranges 
of the peak intensity variations during the soliton collision vs. 
relative velocity shown with shading; dashed line gives the 
maximum possible value. (c,d) Intensity profiles illustrating 
soliton collisions at small (/ii — 0.5) and large (/ii = 5) rela- 
tive velocities. 

Finally, we note that these results can be applied to 
bright solitons existing on top of a background , and 
also to dark solitons which can exist in media with a 
self-defocusing Kerr-type nonlinearity Indeed, it 

has recently been demonstrated [O that the total in- 



tensity profile is uniquely determined by the eigenvalues 
of the fundamental solitons, and does not depend on the 
angular distribution of radiation modes which form the 
background. Therefore, in these cases, our estimates are 
valid for the maximum changes of intensity relative to 
the background level /;,, i.e. maxa, \I{x, z) — Ib\ < hnf- 

In conclusion, we have obtained an accurate estimate 
for the peak intensities of multi-soliton complexes for a 
homogeneous medium with a Kerr-type nonlinearity in 
the (l-l-l)-dimensional situation. We have demonstrated 
that incoherent coupling can result in a decrease of peak 
intensities when the relative soliton velocities are small 
(or zero) and nonlinear interaction is strong. On the 
other hand, solitons with large relative velocities, where 
the interaction is weak, can roughly be superimposed as 
linear waves. 

The authors are members of the Australian Photonics 
Co-operative Research Centre (APCRC). We are grate- 
ful to Adrian Ankiewicz for a critical reading of this 
manuscript. 
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